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ABSTRACT 

A technique for determining flight director laws for the 
longitudinal control of a V/STOL aircraft in landing 
approach is evaluated. The method is based on the application 
of an optimal control model for the human pilot. The vehicle 
studied was the UH-1H helicopter at three approach ground- 
speeds: 60 knots, 40 knots, and 20 knots. The two pilot 
outputs were longitudinal cyclic and collective. In the 
analysis, ten pilot "transfer functions" which relate the 
two control variables to the five displayed and perceived 
quantities were obtained, These transfer functions were 
then utilized to obtain the respective flight director laws. 
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INTRODUCTION 



I . 

VTOL aircraft have received considerable attention for 
both military and civilian use, due to their unique ability 
to operate in areas not accessible to conventional aircraft. 

The rotary-wing aircraft, in particular, has shown such 
increased diversity that display systems which will allow 
VTOL operations in zero/zero conditions are being actively 
sought. Well-designed flight directors are the heart of 
such systems, 

A flight director is a display system which provides 
control commands to the pilot to enable him to complete 
demanding flight tasks with relative ease and precision. 
Recently, Levison /Ref. If has proposed an analytical 
technique for flight director design based upon the optimal 
pilot model pioneered by Kleinman, Baron, and Levison /Ref. 2/. 
A modified form of this technique is utilized in this thesis 
to obtain the longitudinal flight director laws for the UH-1H 
helicopter in three approach conditions in the presence of 
random vertical and horizontal atmospheric turbulence. Five 
vehicle motion quantities, which are normally directly 
displayed or perceived by the pilot , are blended to drive two 
display symbols, the cyclic director and the collective 
director. The design was undertaken at three different 
groundspeeds : 60 knots, 40 knots, and 20 knots. The pilot- 

vehicle system for the optimal modeling procedure is shown 
graphically in Figure 1, 
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FIGURE 1. PILOT - VEHICLE SYSTEM MODEL 



II. METHOD OF ANALYSIS 



A. THE MODELING HYPOTHESIS 

Subject to his inherent limitations, the well-trained, 
well motivated pilot behaves in an optimal manner. The 
pilot's control characteristics can be modeled by the solution 
of an optimal linear control problem and an optimal estimation 
problem with certain "modifications." 



B. MODIFICATIONS FOR PILOT MODELING 

1. A pure time delay is included in each of the pilot's 
control outputs. 

2. Each output neuromuscular system is modeled as a 
first order lag. 

3. Each observed variable is assumed to contain pilot 
induced additive white noise which scales with the 
variance of the observed variable. Also, each 
control output is assumed to contain pilot induced 
additive white noise which scales with the variance 
of the control motion. 

4. If a variable is displayed explicitly, the pilot 
also perceives the first derivative of the variable 
but no higher derivatives. The first derivative is 
also noise contaminated. 

5. The index of performance for the optimization 
procedure is chosen subjectively to mirror what the 
display system designer believes to be the task and 
control objectives perceived by the pilot. 



C. LONGITUDINAL HELICOPTER EQUATIONS OF MOTION 

The longitudinal helicopter equations of motion are 
shown in state variable form on the next page. The 
assumptions used in the derivation of these equations follow 
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1. The vehicle is idealized as a rigid airframe to which 
is attached a rotor. 

2. The rotor is described by the tip path plane whose 
orientation determined propulsive and aerodynamic 
forces and moments. 

3. No rotor degrees of freedom are considered other than 
control inputs which serve to describe instantaneous 
tip path plane orientation. 

4. All coupling between longitudinal and lateral 
motion is ignored. 

5. There is linearized small perturbation motion about 
the horizontal reference flight path. 



D. ATMOSPHERIC TURBULENCE 

The power spectral densities for the atmospheric turbulence 
are those suggested by Hart, Adkins, and Lacau /Ref. 3 /, 
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Using the concept of a white noise excited shaping filter, 
the following state equations can be developed. 
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( 11 ) 



u = -$,u + a,v- 

g 2 g 22 

Equations (10) and (11) now augment equations (1) - (5), 

E, NEUROMUSCULAR AND TIME DELAY EQUATIONS 

As noted in section II. B. , a pure time delay and first 
ordger lag are included in the pilot's output. The delay is 
approximated by a second-order Pade ' function. The quality 
of first and second order Pade' approximations are indicated 
in Figure 2 for the frequency interval of interest in pilot 
modeling, ojt < 5 RAD. For ease of development, a portion of 
Figure 1 is shown here. 
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FIGURE 2. PADE' COMPARISONS 



15 



( 14 ) 



.. 8 ♦ ^ 16 

U 2 d + x U 2 d + T 2 U 2 d = U 2 - 



8 • 16 
T U 2 + 72 u 2 



Now let 



d i = u id ~ u i 



. h ^ 16 
d 2 = d ! + — U 1 



d = d - — u 

1 2 T 1 



(15) 

(16) 
(17) 



d = - 



16 



8 



128 



— 5 - d , - — d, + , 

2 1 ^ 2 ^2 



U , 



a 3 v 3 



t d l T + t ^ l T 

B a n a n b a n a n 



(18) 

(19) 



Using the same technique to develop d , d , and <$ as was 

’ ** C 

used to develop d t , d 2 , and 6 ^, one obtains: 



a 3 = d„ - ^ u 2 



(20) 




Equations (1) - (5), (10) - (11), and (15) - (22) now 
constitiute the state description of the pilot-vehicle 
system. 
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F. DISPLAYED VARIABLES 

The displayed variables listed below are those quantities 
which are assumed to be displayed or perceived by the pilot, 
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- Displayed Groundspeed Deviation = KjU 



P 



P 



1 



2 



3 



4 



5 



- Perceived Time Rate of Change of Glideslope 
Deviation = K 2 (-w + U Q 0 ) 

- Perceived Pitch Rate = K q 

3 

= Displayed Pitch Angle Deviation = K l) 0 

- Displayed Glideslope Deviation = K $ h 



The are display gains. For example: 

„ _ | „ | radians subtended at the pilot's eye by display 

3 ' 3 ■ rad/sec pitch 



element motion 
rate 

Table 1 lists the display gains for this study. 

According to modification (4) in section II, B., u 
should also be a perceived variable. However, since all 
perceived variables will be used in generating the director 
laws to be discussed, they must be measurable. If u repre- 
sents airspeed, u will be difficult to sense. For this 
reason u was not considered a perceived variable (despite 
the fact the u represented groundspeed in this particular 
analysis ) , 

It should be emphasized that neglecting u imposes little 
constraint on the model's validity since this variable 
is associated with the phugoid and will be quite small 
throughout the approach. 
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Display Gain 


Value 


K 1 


3.44 x 10' 4 rad/ ft/sec 


K 2 


6.67 x 10- 4 rad/ ft/se( . 


K 3 


9.55 x 10~ 2 rad/ , , 

' rad/ sec 


K 4 


9.55 x 1CT 2 rad/ , , 

' rad/sec 


K 5 


6.67 x 10“ 4 rad/ f 



TABLE I. Display Gains 

G. OBSERVATION AND MOTOR NOISE 

At present it is not possible to separate experimentally 
the various sources of pilot "remnant." In this paper, as 
is done in reference 2, observation noise and motor noise 
were taken to be representations of remnant. Observation 
noise was taken to be pilot injected noise in observing the 
displayed quantities, and motor noise was considered to be 
random error in control movement execution. Both were 
considered independent, Gaussian, white noise processes. 

The pilot related noise levels were set at values 
considerably larger than those found in documented laboratory 
experiments (e.g. Ref. 2) so that the design would be less 
sensitive to pilot noise and more forgiving of actual non- 
optimal pilot behavior. Consequently, numerical values for 
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the observation noise and motor noise were chosen as p = 0,1, 
and p' = .01. 

H. SYSTEM EQUATIONS AND MATRIX NOTATIONS 

The system equations listed below define the optimal state 
feedback controller and estimator problem. A thorough develop 



ment of these equations can be found in reference 4. 

x(t) = A x (t) + B u (t) + y w (t) (23) 

£(t) = £ x (t) (24) 

z(t) = II x (t) + v (t) = z (t) + v (t) (25) 

P 

E £ w( t )w T (T + i)j = F 6 (t) (26) 

E £v(t)v T (T + x)j = G 6 (t ) (27) 

T 

J = f / [ Z T ( t ) Q y (t) + u T (t) R u (t)J dt (28) 



o 

where 

A is an n x n plant matrix 
x(t) is an n x 1 state vector 
B is an n x p control matrix 
u(t) is a p x 1 control vector 
Y is an n x t disturbance matrix 
w(t) is a t x 1 disturbance vector 
^_(t) is a q x 1 output vector 
c is a q x n output matrix 

Here, w(t) is a vector of linearly uncorrelated, zero 
mean white noise signals with Gaussian amplitude probability 
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distribution functions, The elements of w(t) are assumed to 
be sample functions from n random processes which are each 
ergodic and are jointly ergodic. The covariance matrix for 
w(t) is 

E £ w(t ) w T (t + x )] = F S(t ) (29) 

where 6(x) is the unit impulse function. 

The measured quantities on sensor signals are 

z(t) = H w(t) + v (t) (30) 

where 

z(t) is a u x 1 measurement vector 
H is a u x n measurement matrix 
v(t) is a u x 1 measurement noise vector 
The elements of v(t) are assumed to be sample functions from 
p random processes each of which are ergodic and jointly 
ergodic, The covariance matrix for v(t) is 

E [v(t) v T (t + t)J = G 6(t) (31) 

The system is assumed to be completely controllable and 
completely observable. It is desired to find the control 
function u(t) which minimizes the quadratic scalar index of 
performance 

T 

J = tIZ f J* £ y T (t ) Q y(t) + u T (t) R u(t)J dt (32) 

o 

where 
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Q is a q x q symmetric output cost weighting matrix and 
at least positive semidefinite 



R is a p x p symmetric control cost weighting matrix and 
positive definite 



The solution to the linear quadratic Gaussian control 
problem can be outlined as follows: 

a. ) The optimization problem can, by the called Separation 

Theorem, be broken up into two separate problems, an 
optimal control problem and an optimal estimation 
or filtering problem. 

b. ) The optimal estimation or filtering problem generates 

an optimal estimate, x(t) of the state x(t). This 
estimate is optimal in the sense that 

T 

lim 1 r x^( t ) x(t) dt 

T->a> T / 

o 

is minimized, where x(t) is the estimation error defined as 



x(t) = x(t) - x(t) 



(33) 



The optimal estimator (or Kalman filter) has the form 



&(t) = A x(t) + B u(t) + K |^z(t) - H x(t) 

The estimator gains are given by 

T -1 

K = P H G 



(34) 



(35) 



where P is the error covariance matrix 

E J^x(t)x ,r (t + x)j = P 6(t) (36) 

P is the positive definite solution to the steady - state 
filter matrix Riccat i equation 



A P + P A 



T T -1 

yfy -phg- l hp = o 



(37) 
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c,) The optimal control problem generates an optimal 

control law u(t) which is a linear function of the 
estimated state 



u(t) = - L x(t) (38) 

where L is a p x n optimal controller gain matrix. The gain 
matrix L is identical to the one obtained by solving the 
optimal control problem with no system disturbance, exact 
state information, and the index of performance given by 

J = [y T (t) Q £(t) + u T (t) R U(t)l dt (39) 

the controller gain matrix L is given by 

L = FT 1 B T S (40) 

where S is the positive definite solution to the steady-state 
control matrix Riccat i equation 

-SA-A T S-C T QC+SB R -1 B T S = 0 (41) 

It can be shown that the state covariance matrix 

E |jc(t) x T (t + T) ] = (P + M) 6(t) (42) 

where P is the solution ot the filer matrix Riccati equation 
and M is the positive definite solution to 

(A - B L) M + M (A - B L) T + K G K T = 0 (43) 

In addition to the solutions outlined above, it can be 
shown that the transfer matrix relating the Laplace transform 



22 



of the optimal control law u(t) to the Laplace transform of 
the measurement vector z(t) (with v (t) = 0) is given by 



U (S) = - L (SI_ - A + B L + K H)" 1 K Z (S) 



(44) 



where 

u (S) =J[u (t)] (45) 

Z (S) = /[z (t)J (46) 

The state variables for this study were chosen in the 
following order 



u 

w 

q 

g 

h 




The following matrix tables were developed from equations 
(1) - (5), (10) - (11), and (17) - (22). The matrices are 
labeled in accordance with equations (23) - (27), 

As Table VI indicates, all but one of the elements of 
F and G are dependent upon variances of system variables 
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which are not known a-priori , Estimates of these variances 



must be made, the solution to the optimal estimation and 

control problem obtained, and the resulting variances used 

in a second iteration. This iterative process continues 

until the equations for the a. and V in Table VI are 

i 

satisfied. As reference 2 points out, 2*m iterations are 
usually involved, where m is the number of displayed and 
perceived quantities. 

The output vector y utilized in the index of performance 
is given in Table V as 



y = 



u 

q 

h 

i - 



The elements of the index of performance weighting 
matrices Q and R were chosen as the reciprocals of the 
"maximum allowable" deviations of the output and controls. 
Thus, when an output or control variable attains these 
subjectively chosen magnitudes, it makes a contribution of 
unity to the integrand of the index of performance. The 
maximum allowable deviations were chosen subjectively and 
varied with the flight conditions as Table V indicates. 



24 



A MATRIX 
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TABLE II 
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TABLE III. 
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Matrix C Matrix 20 Knots 
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TABLE VI. 



Ill, PILOT MODELING EXAMPLE 



This section presents a numerical illustration of the 
pilot modeling technique developed in Section H. The 
longitudinal dynamics of a UH-1H helicopter at an approach 
groundspeed of 60 knots are utilized. The normalized UH-1H 
longitudinal derivatives, in a stability axis system, are 
shown in Table VII, 

A modified form of the Variable Automatic Synthesis 
Program (VASP) /Ref. 5/ was utilized to solve the optimal 
estimation and control problem. The solution to this 
problem constitutes the pilot-vehicle model. After roughly 
10 iterations the solution for the 60-knot case converged. 
Table VIII shows the root-mean-square (RMS) performance 
figures. The ten pilot transfer functions relating the cyclic 
and collective control variables to the five displayed and 
perceived quantities were obtained. Figures 3 and 4 are 

6 b 6 c 

Bode plot representations for the — (s) and — (s) pilot 

u h 

transfer functions respectively. 

Pure gain approximations were then made to each of the 
transfer functions in the frequency range of interest for 
modeling purposes: . 1 < w < rad/sec. The gains were then 

normalized by dividing by the magnitude of the largest gain 
for each control, The values in Table IX resulted. For 
example, using the gains of Table IX the cyclic and collective 
flight director laws become 



30 



( 47 ) 



6„(t) = - 8,42 • l(T 4 u(t) .. 2,45 « l(T 4 £i(t) + q(t) 

+ ,1945 Q(t) - 4.02 • 10" 4 h(t) 

6 c (t) = 2,47 * 10“ 4 u(t) - 5.52 ♦ 10~ 3 h(t) - q(t) 

- ,1358 0(t) + 1.96 * 10" 3 h(t) 
where S' and S' are the cyclic and collective director 

D U 

signals respectively. 



(48) 
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The NORMALIZED UH-1H LONGITUDINAL DERIVATIVES IN STABILITY AXIS SYSTEM 



TABLE VII 
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71803E 00 1/ft- 



RMS 




Performance 


Vu 2 (t) 


= 


3,03 ft/sec 




= 


5,11 ft/sec 


V q 2 (t ) 


= 


1,01 degree/sec 


Ve 2 (t) 


= 


1 . 46 degree 


Vh 2 ( t) 


= 


9,32 ft 


g (t) 


= 


5 ft/sec 


Vu 2 (t) 

g 


= 


5 ft/sec 


V6 g 2 (t) 


— 


1,12 in 


Vi c 2 ( t) 


- 


1,35 in 



TABLE VIII . 
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Frequency (rad/sec) 



FIGURE 4, 



TRANSFER FUNCTION — - 

h 
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DIRECTOR 


SENSORY 

VARIABLE 


LOW -FREQUENCY 
GAIN 




Airspeed 


-.000842 sec 




Sink Rate 


-.000245 sec 


CYCLIC 


Pitch Rate 


.. ft-sec 
rad 




Pitch 


,1945 IX , 
rad 




Height 


-.000402 ~ 




Airspeed 


,000247 sec 




Sink Rate 


,00552 sec 


COLLECTIVE 


Pitch Rate 


1 ft-sec 
rad 




Pitch 


ft 

-,1358 IX , 
rad 




Height 


,00196 



TABLE IX, Gains for UH-1H Helicopter Director Laws - 
Velocity 60 Knots 
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IV, RESULTS AND CONCLUSIONS 



The normalized director gains for the 40 and 20-knot 
approach speeds are shown in the following tables. It should 
be noted that all of the Bode magnitude plots for the thirty 
pilot transfer functions (ten each for each of the three 
approach speeds) could be approximated by a fifth-order low 
pass filter. The break frequency varied somewhat between 
0.4 and 0.5 rad/sec. Only the director control gain changed 
from function to function. 

In order to ascertain the reason for the rather dramatic 
sign reversals which occurred in some of the larger gains as 
the flight condition changed (e.g. in the collective-to-pitch 
rate and collective-to-pitch gains of Table IX and Table X), 
the author re-ran the 40-knot case with identical pilot 
parameters, specifically Q and R matrices, as in the 60-knot 
case. The results were similar to those of Table X. This 
indicated that the flight condition dictated the gain sign 
variation and not the subjective selection of the Q and R 
matrices in the pilot model. 

The director laws implicit in Tables IX - XI must be 
evaluated in piloted simulation before the efficiency of the 
design method outlined in this thesis can be determined, 
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DIRECTOR 


SENSORY 

VARIABLE 


LOW-FREQUENCY 

GAIN 




Airspeed 


.000536 sec 




Sink Rate 


-.005287 sec 


CYCLIC 


Pitch Rate 


- „ ft-sec 
,U rad 




Pitch 


.374 ft , 
rad 




Height 


-.00175 ~ 
It 




Airspeed 


.000426 sec 




Sink Rate 


-.0044 sec 


COLLECTIVE 


Pitch Rate 


- _ ft-sec 
,U rad 




Pitch 


,41962 ft , 
rad 




Height 


-.00188 



TABLE X. Gains for UH-1II Helicopter Director Laws - 
Velocity 40 Knots. 
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DIRECTOR 


SENSORY 

VARIABLE 


LOW-FREQUENCY 

GAIN 




Airspeed 


-.0126 sec 




Sink Rate 


,004448 sec 


CYCLIC 


Pitch Rate 


1 ft-sec 
rad 




Pitch 


ft 

.277 , 

rad 




Height 


,00259 || 




Airspeed 


-.01257 sec 




Sink Rate 


.00435 sec 


COLLECTIVE 


Pitch Rate 


^ ft-sec 
rad 




Pitch 


f t 

.3005 ii-r 
rad 




Height 


.002535 || 



TABLE XI , Gains for UH-1H Helicopter Director Laws - 
Velocity 20 Knots. 
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